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Abstract. Grover's database search algorithm, although discovered in the context of quantum 
\ computation, can be implemented using any physical system that allows superposition of states. A 

■ physical realization of this algorithm is described using coupled simple harmonic oscillators, which 
' can be exactly solved in both classical and quantum domains. Classical wave algorithms are far 
, more stable against decoherence compared to their quantum counterparts. In addition to providing 

(N . convenient demonstration models, they may have a role in practical situations, such as catalysis. 

^SJ ' Keywords: Grover's database search algorithm. Catalysis, Resonance 

>■ PACS: 03.67.Lx, 34.10.+x,46.40.-f 

(N : 
^ ■ 

g : COMPUTING WITH WAVES 

o : 

^ ' Any physical system — with some initial state, some final state, and some interaction 

in between — is a candidate for processing information. One only needs to construct a 
suitable map between physical properties of the system and abstract mathematical vari- 
ables. Most of the development in computer algorithms has been in the framework of 

■ "particle-like" discrete digital languages. It is known that "wave-like" analogue compu- 
^ ! tation can also be carried out (e.g. using RLC circuits), but that has not been explored 
qh| as intensively. The obvious reason is that discrete variables allow a degree of precision, 

• by implementation of error correction procedures, that continuous variables cannot pro- 

• T-j . vide. In addition, computational complexity is believed to be the same for digital and 

I analogue algorithms, so the choice between the two is left to considerations of hardware 

(73 ' stability. 

With the advent of quantum computation, several quantum algorithms have been dis- 
covered, which are superior to their counterparts based on Boolean logic. Naturally, 
with both "particle-like" and "wave-like" behaviour at their disposal, quantum algo- 
rithms cannot do any worse than their classical counterparts. It is routinely stated that 
the simple parallelism provided by superposition and interference of quantum states is 
the key ingredient for the superiority of quantum algorithms. Now both superposition 
and interference are generic features of wave dynamics, and it is worthwhile to investi- 
gate the advantages they bring to an algorithm by exploring implementations based on 
"wave-like" behaviour alone. Of course, the classical wave implementations will be less 
efficient than the fully quantum ones, but they are expected to be much more stable (in 
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particular, they have no entanglement and much weaker decoherence) and so may turn 
out to be useful in specific situations. With this motivation, I study in this work, Grover's 
database search algorithm ^ — a straightforward and yet versatile algorithm that allows 
a variety of physical realisations. 

The Optimal Search Algorithm 

Database search is an elementary computational task, whose efficiency is measured 
in terms of the number of queries one has to make to the database in order to find the 
desired item. In the conventional formulation of the problem, the query is a binary oracle 
(i.e. a Yes/No question). For an unsorted database of A'^ items, using classical Boolean 
logic, one requires on the average {Q) = N binary queries to locate the desired item. The 
number of queries is reduced to (Q) = (A^+ 1) /2, if the search process has a memory so 
that an item rejected once is not picked up again for inspection. 

Grover discovered a search algorithm that, using superposition of states, reduces the 
number of required queries to Q = 0{\/N) 111]. This algorithm starts with a superposition 
state, where each item has an equal probability to get picked, and evolves it to a target 
state where only the desired item can get picked. Following Dirac's notation, and using 
the index i to label the items, the starting and the target state satisfy 

\{i\s)\^=\/N, \m^=du . (1) 

The algorithm evolves 1^) towards \t),hy discrete rotations in the two-dimensional space 
formed by I*) and \t), using the two reflection operators, 

Ut = l-2\t){t\, Us=l-2\s){s\ , (2) 

i-UsUtfls) = \t) . (3) 

Ut is the binary oracle which flips the sign of the target state amplitude, while —Us 
performs the reflection-in-the-mean operation. Solution to Eq.© determines the number 
of queries as 

(2e+l)sin-i(l/VA^) = ;r/2 . (4) 

(In practice, Q must be an integer, while Eq.® may not have an integer solution. In such 
cases, the algorithm is stopped when the state has evolved sufficiently close to, although 
not exactly equal to, \t). Then one finds the desired item with a high probability.) 

In the qubit implementation of the algorithm, one chooses N = 2" and the items in the 
database are labeled with binary digits. With the uniform superposition as the starting 
state, 

(/|5) = l/v^, Us = H®"{l-2\0){0\)H®" , (5) 

(H is the Hadamard operator), the implementation requires only 0{log2N) spatial re- 
sources [1]. It has been proved that this is the optimal algorithm for unsorted database 
search [2]. 

Different physical realizations of the database items (e.g. binary labels or individual 
modes) and the target query oracle (e.g discrete binary operation or continuous time evo- 
lution) produce a variety of implementations of this algorithm. In the original version. 



the states are encoded in an n-qubit register, and the oracle is a discrete binary operation 
(denoted by Ut above). In the analogue version of the algorithm, the discrete unitary 
oracle is traded for a continuous time interaction Hamiltonian, which evolves the target 
state somewhat differently than the rest and acts for the entire duration of the algorithm, 
and the number of queries is replaced by the time one has to wait for before finding 
the target state [3]. The wave version of the algorithm requires distinct wave modes, 
instead of n qubits, but does not involve quantum entanglement at any stage [4]. Such a 
wave search has been experimentally implemented using classical Fourier optics, with 
a phase-shift plate providing the oracle [5]. An analogue version of the algorithm has 
been described using a classical coupled pendulum model, where one of the pendulums 
is slightly different than the rest and the uniform superposition state I*) is identified with 
the center-of-mass mode [@]. In what follows, I describe a binary oracle version of the 
wave search algorithm using identical coupled harmonic oscillators. 

HARMONIC OSCILLATOR IMPLEMENTATION 

A harmonic oscillator is the favourite model of physicists. It provides the first approx- 
imation in a wide variety of physical phenomena involving small fluctuations about an 
equilibrium configuration. Its mathematical description contains only quadratic forms, 
in position as well as momentum coordinates. It can be solved exactly in both classical 
and quantum domains, which makes it extremely useful in situations where a cross-over 
between classical and quantum behaviour is to be analyzed. We shall first look at the 
classical system, and then observe that the quantum system essentially follows the same 
pattern. 



Classical Oscillators 

Let the items in the database be represented by identical harmonic oscillators. 
While they are oscillating in a specific manner, someone taps one of the oscillators (i.e. 
elastically reflects it by touching it). The task is to identify which of the oscillators has 
been tapped, without looking at the tapping — quite like a magician who can tell which 
one of his cards was touched when he was blindfolded. The optimization criterion is to 
design the system of oscillators, and their initial state, so as to make the identification as 
quickly as possible. 

Grover's algorithm requires identical coupling between any pair of oscillators. This 
can be accomplished by coupling all the oscillators to a big oscillator, as shown in 
Fig.l. The big oscillator is thus coupled to the centre-of-mass mode, and becomes an 
intermediary between any pair of oscillators with the same strength. The Lagrangian for 
the whole system is 

1 1 ^1 1 

^ = -MX^ - -KX^ + £ - -k{x, -Xf]. (6) 



With the center-of-mass displacement, x = T!iL\^i/N, the Lagrangian can be rewritten 
as 

^ = ^MX^-Ux^ + ^Nmi^-^Nk{x-Xf + f^[^^^ . (7) 

Now we can fix the oscillator parameters to implement Grover's algorithm. In the 
algorithm, we are interested in the dynamics of the tapped oscillator. All the other 
oscillators (i.e. i 7^ influence the dynamics of Xt only through the combination x. The 
dynamics of (A'^ — 2) linearly independent modes orthogonal to Xr and x (they all have 
the form {xj^j —Xj^^t)) decouples from the modes of interest; we can drop them and 
effectively work in the 3-dimensional space of the modes {X,x,Xt}. (In what follows, I 
shall first specify initial conditions such that all Xj^j are identical and all the decoupled 
modes vanish. Subsequently, we will look at the general situation by adding back all the 
decoupled modes.) 

Choosing units of mass and time such that m = \ ,k= 1 , and in terms of the variables 
Y = VmX , y = VNx ^ yt=xt-x, (8) 
the effective Lagrangian becomes 

The potential energy terms in ^eff are easily diagonalized, and yield the eigenvalues 

^2V M ) \ A\ M ) M' 

9 9 K -\- TV 99 K 
(o^ + (oi = l-\ , (0^(01 = — , ft)( = 1 . (10) 
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FIGURE 1. A system of identical harmonic oscillators, coupled to a big oscillator via the center-of- 
mass mode. 



The corresponding eigenmodes are 



{\-col)Y+^y={\-(ol)^/MX + ^x, (11) 



e± - 
et = yt=xt-x . 

The initial uniform superposition state can be realized as all the oscillators moving 
together, while the big oscillator is at rest. 

^ = 0: X=0, X = 0, JC/ = 0, i,=A. (12) 

(We will consider situations with general initial conditions later.) The reflection op- 
erators correspond to shifting the appropriate oscillator phases by half a period. The 
binary tapping oracle can be realized as the elastic reflection illustrated in Fig. 2. That 
implements Ut in the velocity space, by reversing the target oscillator velocity at the 
instance when all the displacements vanish. Time evolution of the coupled oscillators 
redistributes the total kinetic energy, and that can implement the operator Us with a suit- 
able choice of time interval and frequencies. 

With the natural frequency of individual oscillators CO = ^Jk/m = 1, the reflection-in- 
the-mean operation requires 0}± to be rational numbers. Optimization means that they 
should be selected to make the dynamics of the whole system of oscillators have as small 
a period as possible. The solution is not unique. One set of solutions is {p is a positive 
integer) 

2p + \ 1 \6Nm {2p + \fNk 

^+ = ^' ^- = 2 ^ ^=3(2;. + 3)(2;.-l)' ^=3i2iT3)i2i^)- ^''^ 

In these cases, in absence of oracles, the dynamics of the whole system of oscillators 
has the period, T = Ak. The big oscillator returns to its initial rest state (X = 0,X = 0), 
whenever t is an integral multiple of lit, i.e. after every half a period. Time evolution 
for the same interval of half a period reverses x, while leaving Xt ~x unchanged, i.e. it 
implements the operator Us in the velocity space. Thus Grover's algorithm, Eq.©, can 
be realized by applying the tapping oracle at every time interval A? = 2;r. 
A more interesting set of solutions is 

Nm 

co+ = 2p, (0-=0 ^ M= -, K = 0. (14) 

(2p + l)(2p- 1) 



(a) 'ooooirQ) Xt = 0, i:^ > : Just before the oracle 
(6) 'TOsojrQ^ xt = 0, Xt ^ —xt : Elastic reflection 
(c) ^mmi^ Xt = 0, Xt < : Just after the oracle 



FIGURE 2. The binary tapping oracle flips flie sign of the target osciUator velocity, when its displace- 
ment is zero. 



Under these conditions, the big oscillator is not coupled to any support, and e becomes 
a translation mode for the whole system of oscillators. The translation mode can be 
eliminated from the dynamics with the initial conditions 



t = 0: X = 0, X = -—A , Xi = 0, Xi=A. 



(15) 



Then, in absence of oracles, the dynamics of the whole system of oscillators has the 
smallest possible period, T = 2n. After half a period, the big oscillator is back to 
its initial state, x also returns to its initial value, while if — I changes its sign. This 
is equivalent to applying —Us in the velocity space, and Grover's algorithm can be 
implemented by tapping the target oscillator at every time interval A? = Ti. 

There is an important physical distinction between the quantum and the wave im- 
plementations of the amplitude amplification process in Grover's algorithm — quantum 
probability is mapped to wave energy. The enhancement of the quantum amplitude in- 
creases the probability of finding the target state A^-fold, while the enhancement of the 
wave amplitude increases the energy of the target oscillator AT-fold. The well-known 
phenomenon of "beats" is responsible for energy transfer amongst coupled oscillators. 
The elastic reflection oracle does not change energy, and it is interesting to observe that 
the oscillator which is obstructed by tapping picks up energy. 



Now we can look at the behaviour of the wave implementation under more general 
circumstances. First consider the initial conditions. Despite appearances, precise syn- 
chronization of oscillators is not an issue in the algorithm, because of the explicit cou- 
pling to the center-of-mass mode. For instance, the algorithm can be started off with 
an initial push to the big oscillator, X = B,Xi = 0, and the system of oscillators would 
evolve to the stage X = 0,Xi = A. Furthermore, any arbitrary distribution of initial ve- 
locities of oscillators can be accommodated in the analysis by bringing back the (A^ — 2) 
decoupled modes. The decoupled modes have no effect whatsoever on the dynamics of 
the {X,x,Xt} modes. Consequently, they modify the algorithm only to the extent that the 
energy amplification of the target oscillator is limited to the initial energy present in the 
{X,x,Xt} modes, instead of being A^-fold. Explicitly, the maximum gain is 



which can be substantial for the generic situation where the initial Xt and x are of the 
same order of magnitude. 

To extract the maximum gain, the algorithm must be stopped at a precise instant 
(i.e after a precise number of tapping oracles Q); otherwise the evolution continues in 
repetitive cycles. The state evolution in Grover's algorithm is a rotation at a uniform rate 
in the two dimensional \s)-\t) subspace. The average overlap of the target state, with the 
state 1^(2)) after Q queries, is therefore 



Stability Considerations 




(16) 



|(?(e)k)lav=(sin'0)av=l/2. 



(17) 



Thus if the algorithm is stopped at a random instant, the energy gain on the average is 
half of its maximum value in Eq.(fT6l) — which can still be substantially larger than 1 . 

Next consider the effect of damping. Damping for a harmonic oscillator has to be 
analyzed using its equation of motion; it cannot be described using a time-independent 
Lagrangian. The standard description is: 



jc + 27i+£t)oJc = ^ x = Ae ^ cos y^J co^ - t + (j) j . (18) 

The crucial ingredient in the algorithm is the coherence amongst the phases of the os- 
cillators. That is governed by the frequencies of the oscillators, and is independent of 
the amplitudes. For a weakly damped oscillator, its amplitude changes linearly with the 
damping coefficient, while its frequency changes quadratically. The time evolution of 
the above implementation, therefore, remains essentially unaffected if the oscillators ex- 
perience a small damping. The dominant effect is a decrease in the energy amplification 
due to decaying amplitudes. 

Among other variations, simultaneous scaling of masses and spring constants of the 
oscillators (i.e. m,- = tt/m and kj = Uik) does not alter the algorithm at all, since the scale 
factors can be absorbed by redefining x/. One can also contemplate interchange of the 
initial and the final states, since the algorithm is fully reversible. The evolution would 
then run backwards, and the physical interpretation would be the dissipation of a large 
initial energy of the target oscillator in to a uniform distribution over all the coupled 
oscillators. 



Quantum Domain 

The dynamics of harmonic oscillators is simple enough to permit exact quantum 
analysis as well. It is convenient to interpolate between classical and quantum domains 
using the coherent state formulation (see for example, Ref. [7]). Coherent states are 
superpositions of the energy eigenstates, parametrized by a single complex variable a, 

|a) = e-l«lV2y_^|„). (19) 
« \/n\ 



They describe Gaussian wavepackets with minimal spread (i.e. displaced versions of the 
ground state eigenfunction). 



, h . mh(0 

A coherent state with the initial condition a{t = 0) = (Zo h^s energy h(o{\(Xo\^ + j), and 
the centre of its wavepacket performs the same simple harmonic motion as a classical 
particle would. 
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FIGURE 3. Evolution of the coherent state wavefunction of the tapped oscillator, with the initial 
condition ocq = —a. The left half of the figure shows the actual wavepacket in the half-harmonic oscillator 
potential, with the impenetrable wall at x = 0. The right half of the figure shows the image wavepacket 
that ensures the node of the wavefunction at x = 0. For t = nco /2, the wavefunction is purely imaginary, 
but the factor of i is omitted for convenience in drawing the figure. The wavepacket at / = n(0 includes 
the geometric phase of — 1 arising from reflection. 



The wavefunction of the state evolves according to 

|V^(0)) = |ao) =^ \xif{t)) = e~''''/^\aoe-'">'), (22) 
while the explicit structure of the wavepacket is given by 




Thus the classical analysis of previous sections can be carried over unchanged to 
the quantum domain, provided we can figure out how the tapping oracle works for 
the coherent states. The tapped oscillator corresponds to a particle moving in the half- 



harmonic oscillator potential 

V{x) = \kx^ for jc < , y (x) = oo for X > . (24) 

The impenetrable wall at x = is equivalent to enforcing the boundary condition ^/(a: = 
0) = 0. So the energy eigenstates of the half-harmonic oscillator are the same as those 
for the harmonic oscillator, with odd n. It is straightforward to ensure the node at x = 
using the method of images, and the tapped oscillator coherent states become 

|a,)=C(|a)-|-a)) , (25) 

with a time-independent normalization constantC= (1 — e^^l"!')^^/^. Tapping amounts 
to change-over between |a) and \ — a), which reverses {x) and {p) compared to the 
untapped motion. In addition, the wavefunction changes sign, which is the geometric 
phase corresponding to wave reflection. The evolution of a coherent state wavepacket 
undergoing a reflection from the wall is depicted in Fig. 3. 

POSSIBLE APPLICATIONS 

The oscillator based search process discussed above has the same algorithmic efficiency 
as the proposals of Refs.[@, [si 0] — where it differs from them is in the actual physical 
implementation. All these wave implementations require exponentially more spatial 
resources compared to their digital counterparts, 0{N) vs. 0(log2 A'^). On the other hand, 
they reduce the number of oracle calls by exploiting superposition of states. Note that 
no algorithm based on Boolean logic, either with serial or with parallel implementation, 
can reduce the number of oracle calls for an unsorted database search to less than 0{N) . 

Quantum algorithms are superior to wave algorithms, because they can use super- 
position as well as reduce spatial resources. The reduction of spatial resources, how- 
ever, comes with the cost that quantum algorithms have to work with entangled states. 
Quantum entanglement is far more sensitive to decoherence caused by environmental 
disturbances than mere superposition, and that has made physical implementations of 
quantum algorithms very difficult. On the other hand, superposition of classical waves 
can be fairly stable, even in presence of a small amount of damping, and that can make 
wave implementations advantageous in specific physical contexts. 

These features indicate that wave algorithms fall in a regime in between classical 
and quantum algorithms — more efficient than the former and more robust than the 
latter. They are likely to be useful in practical situations, where is not very large and 
environmental disturbances are not negligible. Indeed it is worthwhile to systematically 
explore them, just like randomized algorithms have been [8]. 

In the specific case of the unsorted database search problem, the remarkable simplicity 
of the oscillator implementation would allow construction of convenient demonstration 
models, which can demonstrate the power of quantum algorithms even at school level. 
Going further, the practically useful property of the wave search algorithm is that it fo- 
cuses energy of many oscillator modes into one of them. So it is really interesting to think 
of situations where the energy amplification process can have dramatic consequences. 
One such possibility is very familiar to all of us — the Boltzmann factor where the energy 



is in the exponent — and I describe below a scenario for catalysis where involvement of 
new mechanisms can enhance our understanding of the observed phenomena. 



Catalysis 

There exist a large number of chemical reactions which, although not forbidden by 
energy conservation, are extremely slow because they have to pass through an intermedi- 
ate state of high energy. In these reactions, the dominant term governing the reaction rate 
is the Boltzmann factor, exp{—Ei,/kT), with the barrier energy in the exponent. Only a 
tiny fraction of the molecules in the tail of the thermal distribution are energetic enough 
to go over the barrier and complete the reaction. It is known that the rates of many such 
reactions can be enhanced by orders of magnitude by adding suitable catalysts (enzymes 
in case of biochemical reactions) to the reactants. The conventional explanation for the 
reaction rate enhancement, called transition state theory, is that the catalysts lower the 
energy of the intervening barrier by altering the chemical environment of the reactants. 

The preceding analysis of the wave search algorithm suggests another mechanism 
for catalysis. Vibrations and rotations of molecules are ubiquitous harmonic oscillator 
modes. The catalyst can act as the big oscillator and transfer energy of many modes to 
the reactant which faces the energy barrier. For example, the catalytic substrate can 
have many identical molecules of one reactant stuck to it and vibrating, the second 
reactant then comes drifting along and interacts with one of the stuck molecules, that 
molecule picks up energy from its neighbours and the reaction gets completed. The rate 
enhancement results not from lowering of the energy barrier but from increase of the 
reactant energy. In such a scenario, for maximum efficiency, the physical parameters 
(masses and spring constants) need to have specific values. But even without perfectly 
tuned parameters, there can be partial energy focusing that provides useful increase in 
the reaction rate. Whether this mechanism exists among the known catalysts, or whether 
we can design new type of catalysts that use this mechanism, is open to investigation. 

The contributions of the "transition state" and the "energy transfer" mechanisms to a 
catalytic effect are not mutually exclusive. So it is desirable to identify characteristics 
that can tell one of them apart from another. I point out two observable features that can 
distinguish the role of chemical environment from physical waves: 

(1) Isotopic substitution in the reactants changes physical parameters without altering 
chemical properties. The electronic potential is essentially independent of the nuclear 
mass, and mass dependence enters the conventional transition state theory only through 
diffusion and tunneling effects. This mass dependence is rather weak and monotonic. 
On the other hand, the nuclear mass strongly affects vibration and rotation frequen- 
cies, which can substantially alter the energy transfer amongst coupled oscillator modes. 
Also, resonant energy transfer is not monotonic, i.e. it decreases on either side of the op- 
timal parameter value. 

(2) In the transition state theory, the reaction takes place between individual reactant 
molecules. On the other hand, energy transfer is a cooperative phenomenon that cannot 
occur without participation of nearby oscillator modes or molecules with similar prop- 
erties. It would therefore be enhanced by existence of non-reacting but similar chemical 



bonds close to the reaction site, and also by increasing concentration of the reactants. 



Kinetic Isotope Effect 

In the context described above, isotope dependence of reaction rates is a signal of 
involvement of physical (in contrast to chemical) features in the catalytic process. Many 
examples of isotopic dependence of catalytic reaction rates have been discovered, and 
the effect is referred to as the "Westheimer effect" or the "kinetic isotope effect" The 
effect is the largest for substitution of hydrogen by deuterium or tritium, and has been 
extensively studied for the rupture of C-H/C-D/C-T bonds. The conventional transition 
state theory has been found inadequate for a theoretical understanding in several cases 
involving enzymes, and vibrationally enhanced quantum tunneling has been invoked as 
an alternative [10]. 

Specifically, the observed kinetic isotope effects are divided in to two categories. 
The "primary effect" results from isotopic substitution of the hydrogen atom that is ex- 
changed between the reactants. The "secondary effect" arises from isotopic substitution 
of a hydrogen atom that is not exchanged during the reaction but is adjacent to the ex- 
changed one. Empirical models of reaction coordinate transition have been constructed, 
but they require effects of quantum tunneling under the barrier as well as contribution 
from vibrational bond energy to explain the size of the observed effects. In particular, the 
secondary effect cannot be explained at all without invoking coupled dynamics of atoms 
ill]. All these additional contributions enhance the reaction rates beyond their values 
in the classical transition state theory, and hence would be favoured in useful biochem- 
ical processes by natural selection during evolution. The observed results, however, are 
quoted as ratios of reaction rates to eliminate unknown normalization constants. 

Against this backdrop, the coupled oscillator based energy transfer mechanism de- 
scribed in earlier sections offers several novelties: 

(i) The whole analysis is based on first principles. For coupled oscillators, just a single 
reflection can amplify the vibrational energy by a factor of (3 — jj)^, which can be as 
large as 9. Of course, the maximum amplification possible (with more reflections) is A'^. 

(ii) The vibrational modes contributing to catalysis must be soft, i.e. hco = 0{kT), so 
that they can get thermally excited and participate in the dynamics. This will produce a 
characteristic temperature dependence, with the vibrational contribution dropping out at 
low enough temperatures. 

(iii) Resonant energy transfer requires good frequency matching between coupled oscil- 
lator modes, as in Eqs. (ll3ll4|) . Some of that can be inferred from the molecular struc- 
tures (e.g. similar C-H bonds involved in secondary kinetic isotope effect), and more can 
be tested by spectroscopic methods. 

(iv) The wavefunction sign-flip caused by reflection can switch between bonding and 
anti-bonding molecular orbitals (see Fig. 4), and thus help in transfer of atoms. This fea- 
ture is not related to either the mass or the temperature. 

Clearly, the coupled oscillator inspired catalytic mechanism needs to be explored fur- 
ther, with careful modeling of specific reactions. It would be an impressive achievement 
indeed, if the detailed understanding can be used in design of new types of catalysts. 



(a) 




(b) 



FIGURE 4. Formation of molecular bonds by overlap of electron clouds: (a) a binding orbital is formed 
when the two wavefunctions are in phase, (b) an anti-binding orbital results when the two wavefunctions 
have opposite phases. 
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